Contact mechanics: relation between interfacial separation and load 
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I study the contact between a rigid solid with a randomly rough surface and an elastic block with 
a flat surface. I derive a relation between the (average) interfacial separation u and the applied 
normal squeezing pressure p. I show that for non-adhesive interaction and small applied pressure, 
p ~ exp(— w/tto), in good agreement with recent experimental observation. 
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When two elastic solids with rough surfaces are 
squeezed together, the solids will in general not make con- 
tact everywhere in the apparent contact area, but only 
at a distribution of asperity contact spots[l|,[3,y,0]- The 
separation m(x) between the surfaces will vary in a nearly 
random way with the lateral coordinates x = (x, y) in 
the apparent contact area. When the applied squeez- 
ing pressure increases, the average surface separation 
u — (u(x)) will decrease, but in most situations it is 
not possible to squeeze the solids into perfect contact 
corresponding to m = 0. The space between two solids 
has a tremendous influence on many important processes, 
e.g., heat transfer^, contact resistivity 6] , lubrication 7], 
sealing optical interference ... . In this paper I will 
present a very simple theory for the (average) separa- 
tion u as a function of the squeezing pressure p. I will 
show that for randomly rough surfaces at low squeezing 
pressures p ~ exp(— m/wq) where the reference length uq 
depends on the nature of the surface roughness but is in- 
dependent of p, in good agreement with experiments!^. 

We consider the frictionless contact between elastic 
solids with randomly rough surfaces. If z = /ii(x) and 
/i2(x) describe the surface profiles, Ei and E2 are the 
Young's clastic moduli of the two solids and vi and 1^2 
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FIG. 1: An elastic block squeezed against a rigid rough sub- 
strate. The separation between the average plane of the sub- 
strate and the average plane of the lower surface of the block 
is denoted by u. Elastic energy is stored in the block in the 
vicinity of the asperity contact regions. 



the corresponding Poisson ratios, then the elastic con- 
tact problem is equivalent to the contact between a 
rigid solid (substrate) with the roughness profile h(x) = 
/ii(x) -t- /i2(x), in contact with an elastic solid (block) 
with a flat surface and with an Young's modulus E and 
Poisson ratio ly choosen so that [13] 
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Introduce a coordinate system xyz with the xy-plane 
in the average surface plane of the rough substrate, and 
the z-axis pointing away from the substrate, see Fig. [1] 
The separation between the average surface plane of the 
block and the average surface plane of the substrate is de- 
noted by u with u > 0. When the applied squeezing force 
p increases, the separation between the surfaces at the in- 
terface will decrease, and we can consider p = p{u) as a 
function of it. The elastic energy Ua\{u) stored in the sub- 
strate asperity-elastic block contact regions must equal 
to the work done by the external pressure p in displacing 
the lower surface of the block towards the substrate, i.e., 
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/ du'Aop{u') = U,x{u) (1) 
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or 



p{u) 



1 dU,x 
Aq du ' 



(2) 



where is the nominal contact area. Equation (2) is 
exact. Theory shows that for low squeezing pressure, the 
area of real contact A varies linearly with the squeezing 
force pAo, and that the interfacial stress distribution, and 
the size-distribution of contact spots, are independent of 
the squeezing pressure[ll, 12|. That is, with increasing p 
existing contact areas grow and new contact areas form 
in such a way that in the thermodynamic limit (infinite- 
sized system) the quantities referred to above remain un- 
changed. It follows immediately that for small load the 
elastic energy stored in the asperity contact region will 
increase linearly with the load, i.e., Ue\{u) = uoAop{u), 
where uq is a characteristic length which depends on the 
surface roughness (see below) but is independent of the 
squeezing pressure p. Thus, for small pressures (2) takes 
the form 



p{u) = -Mo 
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FIG. 2: The parameters a and /3 as a function of the Hurst 
exponent H for three different values of the ratio qi/qo- 
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in good agreement with experimental data for the contact 
between elastic solids when the adhesional interaction be- 
tween the solids can be neglected [9||. We note that the 
result (3) differs drastically from the prediction of the 
Bush et al theory 11 3| . and the theory of Greenwood and 
Williamson (GW) 14*, which for low squeezing pressures 
(for randomly rough surfaces with Gaussian height dis- 
tribution) predict p{u) ^ u~°exp(— 6u^), where a = 1 in 
the Bush et al theory and a = 5/2 in the GW theory. 
Thus these theories do not correctly describe the interfa- 
cial spacing between contacting solids. 



The elastic energy Ud has been studied in Ref. [15 1 
and I4I, and in the simplest approximation it takes the 
form 



f/ol 
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where E* = E/{1 - v'^) and where P{q) = A{C)/Aq is 
the relative contact area when the interface is studied at 
the magnification C = q/qo- The surface roughness power 
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C{q) = 



d'^x{h{-Kh{0))e' 
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where (..) stands for ensemble average. The parameter 
7 = 1 in the simplest case but in general one expect 7 < 1 
(but of order unity) to take into account that the elastic 
energy stored in the contact region (per unit surface area) 
is less than the average elastic energy (per unit surface 
area) for perfect contact, see Ref. [1^. We will use the 
contact mechanics theory of Persson, where for elastic 
non-adhesive contact the function 17, ISl 
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where s{q) = w{q)/E* with 



dx e 



w{q) = ( TT / dq' q'^C'iq') 
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Using (5) gives 
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Substituting (4) and (7) in (2) gives 
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p{u) = -V^7 / dq g2c(9)w((7)e-['^'(«)P/^*l^ 
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Integrating this from u = (complete contact, corre- 
sponding to p — 00) to u gives 
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For very small pressures we get from (8): 

U = -Uolog(p/pc) 



(8) 



or 



where 



P = Pee 



-u/uo 



uo = V^l / dq q^C{q)w[q) 
Jqa 

and where the cut-off Pc is determined by 
Pc = E*exp{-{\ogw)) 



(9) 
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with 
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If we assume that the substrate surface roughness is self 
affine fractal for qo < q < qi we get[l6j 



C{q)^ 



-2(H+1) 



(13) 



where H = 3 — Df, where Df is the fractal dimension. 
The mean of the square of the substrate surface height 
profile is {h^) = /in„s- Substituting (13) in (10) gives 



where 



1/2 rqi/qo 

7 / dx g{x) 



(14) 



(15) 



with 

g{x) = {x- 

Substituting (13) in (12) gives 

(logw) = -log(go^rms) ^ log/3 

where 

dx g{x)\oA^^ (x2(i-^)-l) 



log/3 



2S1^''^^^dxg{x) 
Substituting (14), (16) and (11) in (9) gives 

or 

u = a^^/ii-mslog iPqoK-^sE* /p) 



(16) 
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In Fig. [5] we show the parameters a (for 7=1) and /3 as a 
function of the Hurst exponent H for three different val- 
ues of the ratio qi/qo- Most surfaces which are self affine 
fractal have the Hurst exponent H > 0.5 (or the fractal 
dimension D{ < 2.5). For such surfaces the parameters a 
and P are nearly independent on the ratio qi /q^ between 
the highest qi and smallest go wavevector included in the 
analysis. Surfaces prepared by crack propagation or by 
bombardment with particles typically have H « 0.8 (or 
the fractal dimension _Dt = 3 — -ffw2.2) for which case 
a ~ I and /? w 0.5. In this case we have 



p w 0.5qoh,^sE*e' 



or 



7/li.mslog(0.5go/irms-E*/_p) 
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FIG. 3: The relation between the (natural) logarithm of the 
squeezing pressure p and the interfacial separation it for an 
elastic solid squeezed against a rigid, self affine fractal sur- 
face with the Hurst exponent H = 0.8. The surface has the 
rms roughness /irms = 6 fim and the upper and lower cutoff 
wavevectors are qi — 7.8 x 10^ and qo — 1 x 10* m~^. 



where we have reintroduced the factor of 7. 

In Fig. [3] we show the relation (for 7 = 1) between 
the (natural) logarithm of the squeezing pressure p and 
the interfacial separation u for an elastic solid squeezed 
against a rigid, self affine fractal surface with the Hurst 
exponent H — 0.8, as calculated directly from (8) by 
numerical integration. The dashed line shows the large- 
distance asymptotic behavior given by (21). Note that 
deviation from this logarithmic relation only occurs for 
u < /irms, and that there is a sharp increase in the squeez- 
ing pressure for u < h^ms- Both these facts are in accor- 
dance with the experimental data presented in Ref. 
Note that the slope of the dashed line in Fig. [3] is close to 
— 1, as expected from (21), and that the dashed curve for 
u = gives log(p/£'*) « -3.8 or p/E* « 0.02, which is 
similar to the prediction of (21): O.Sgo'irms ~ 0.03. I have 
performed calculations of the p{u) relation using several 
measured surface roughness power spectra C{q), but the 
general form is always as in Fig. [3l and the slope of the 
line (using 7 = 1) in the linear region is always close to 
-1. 

Pei et al[l^ have performed a finite element method 
computer simulation of the contact mechanics for a poly- 
mer surface, using the measured surface topography '9'| as 
input, squeezed against a flat surface. They found that 
for large separation p ~ exp(— u/7/ii.ins) with 7 w 0.6, 
which is consistent with our numerical and analytical re- 
sults. 

The theory presented above can be easily generalized 
in various ways. Thus, it is possible to include the ad- 
hesional interaction [iSIIil. In this case the work done 
by the external pressure p will be the sum of the stored 
(asperity induced) elastic energy plus the (negative) ad- 
hesional energy, i.e., the right hand side of (1) will now 



be Uei + Ua,d- The theory can also be apphed to study 
how the spacing u{() depends on the magnification. Here 
is the (average) spacing between the sohds in the 
apparent contact areas observed at the magnification C. 
The quantity u{Q is of crucial importance for lubricated 
seals [2^ • The results of these generalizations of the the- 
ory will be presented elsewhere. 

Finally we note that the observation of an effective 
exponential repulsion has important implications for tri- 
bology, colloid science, powder technology, and materials 
scienceiOj. For example, the density or volume of gran- 
ular materials has long been known to have a logarith- 
mic dependence on the externally applied isotropic pres- 
sure or stress, as found, for example, in the compression 
stage during processing of ceramic materials [isj]. Recent 
work on the confinement of nanoparticles have also in- 
dicated an exponential force upon compression |24l|. sug- 
gesting that this relationship could be prevalent among 
quite different types of heterogeneous surfaces. 
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